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Time-periodic perturbations due to classical electromagnetic fields are useful to engineer the topo¬ 
logical properties of matter using the Floquet theory. Here we investigate the effect of quantized 
electromagnetic fields by focusing on the quantized light-matter interaction on the edge state of a 
quantum spin-Hall insulator. A Dicke-type superradiant phase transition occurs at arbitrary weak 
coupling, the electronic spectrum acquires a finite gap and the resulting ground state manifold is 
topological with Chern number ±1. When the total number of excitations is conserved, a photocur¬ 
rent is generated along the edge, being pseudo-quantized as culn(l/tj) in the low frequency limit, 
and decaying as l/uo for high frequencies with uj the photon frequency. The photon spectral function 
exhibits a clean Goldstone mode, a Higgs like collective mode at the optical gap and the polariton 
continuum. 
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Introduction. The conventional band theory of solids 
has attracted renewed attention after the theoretical pre¬ 
diction and experimental discovery of topological insu¬ 
lators (TIs) [ij-lsl . The topological protection of their 
edge/surface states, together with the spin-momentum 
locking due to the strong spin orbit coupling, yields a va¬ 
riety of novel phenomena, such as magnetoelectric effect, 
axion electrodynamics, surface Hall effect, not to mention 
topological superconductivity and Majorana fermions. In 
addition to their theoretical appeal, they also hold the 
promise to realizing practical quantum computers and 
spintronical devices. 

In spite of the growing interest, materials with non¬ 
trivial topological properties are still scarce. Achieving a 
topological band structure often requires engineering the 
Bloch band structure. Recently, time periodic pertur¬ 
bations, mostly due to electromagnetic fields, have been 
proposed and used[3-S| to manipulate the band structure 
using the Floquet theory, the temporal analogue of Bloch 
states. Although the steady state of the resulting Floquet 
topological insulators often possesses a topology different 
from that of their static parents, the actual occupation 
of these states is still far from being understood^ albeit 
it is an essential ingredient for determining the physi¬ 
cal response of the system. The occupation depends on 
the sources of relaxation, e.g. coupling to heat baths, 
momentum scattering, interaction etc., and the conven¬ 
tional ground-state picture with the minimal energy con¬ 
figuration cannot be used for these driven models. To cir¬ 
cumvent this problem, a classical driving field can be re¬ 
placed by its quantum counterparts having its own quan¬ 
tum dynamics. This idea has been pursued for the cre¬ 
ation of Majorana fermions using not only classical driven 
systems [l^ but rather their quantum counterparts using 
high-Q electromagnetic cavity [lH. 

The properties of classical and quantized electromag¬ 
netic fields have long been investigated. In connec¬ 


tion with their interaction with a single atom 12 , [H, 
e.g., many similarities have been demonstrated between 
the corresponding Rabi and Jaynes-Cummings models. 
Besides the similarities, quantum fields also provide 
novel effects and collective phenomena such as sponta¬ 
neous symmetry breaking and the superradiant phase 
H^ lSl due to light-matter interaction in 


transition 


the Dicke model 16 


iH. 


In this work, we combine TIs and quantized electro¬ 
magnetic fields to explore the quantum mechanical ver¬ 
sion of the model of Ref. 18|, shown to display a Flo¬ 
quet topological phase transition. We demonstrate that 
a superradiant quantum phase transition exists in the 
quantum context, and is accompanied by a transition to 
a topologically non-trivial ground state manifold and by 
an induced photocurrent. 



FIG. 1. (Golor online) Left: a single oj frequency boson mode 
is coupled to a continuum of fermionic degrees of freedom, 
which is a quantum spin-Hall insulator in the present case. 
This is the essential difference between our setup and con¬ 
ventional Dicke/Tavis-Gummings models, where the fermions 
possess a discrete spectrum. The new fermionic spectrum be¬ 
comes gapped and shifted (dashed line). Right: cartoon of 
the experimental setup corresponding to Eq. ([T]) of SHI with 
spin filtered edge states placed into an optical cavity. 


The model. We consider the edge dynamics of a spin- 
Hall insulator [l^ (SHI), interacting with circularly polar¬ 
ized quantized electromagnetic field through the Zeeman 
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term (see Fig. [ij. The Hamiltonian is 

H = uja+a + v'^pS^ + ^'^{a+Sp +aS+) , (1) 

Here v is the Fermi velocity of the edge state, lo and a 
denote the energy and annihilation operator of the cavity 
mode, and Cp,o- annihilates an electron with momentum p 
and spin cr. The last term describes the Zeeman interac¬ 
tion with the spin operators expressed as 5”+ = _|_Cp__ 

and Sp = J2cr=±'^^p a-'^p,<y The Zeeman coupling de¬ 
pends on the photon’s frequency [l^ as 5 = ^Jg\uj\, and 
L is the dimensionless length of the edge [2^ . 

Let us first highlight the connection of Eq. o to 
Floquet topological insulators using the one body pic¬ 
ture of the polaritons. Considering just a single mo¬ 
mentum p, the Hamiltonian reduces to the Jaynes- 
Cummings modelfl^. The spectrum consist of energy 

pairs e±(p, m) = ui{m — i) ± {vp — with m 

bein g a positive integer, plus a single mode e_(p, 0 ) = 
—vp [2l|. These polariton energies are closely related to 
the Floquet spectrum of the same problem in the pres¬ 
ence of classical electromagnetic held Q, when the pho¬ 
ton state is prepared in a coherent state with m 3> 1 . 
Importantly, similar to the classical case, the spectrum 
becomes gapped due to interaction with the quantized 
light held for all m > 0 ( 2 ^, and for each hxed m > 0 , the 
contribution of all e_(p, m) states amounts in an electric 
current identical to that in the classical case, quantized 
as j = in agreement with Floquet adiabatic charge 
pumping. 

This one body description has limited meaning, since 
all spins interact simultaneously with the same bosonic 
mode, and the many-body spectrum involves all of them. 
The Hamiltonian Eq. o is, however, integrablejl^ for 
any number of spins, and coincides with the inhomoge¬ 
neous Dicke model [2^. In addition to the energy, the 
total number of excitations, Q = a+a -I- | Sp is a 
conserved quantity, which generates a 17(1) symmetry. 
Without g, the photon state is unpopulated and the to¬ 
tal spin is unpolarized, therefore (Q) = 0. Note that this 
applies to positive helicity (A) of light, A = 1, for neg¬ 
ative helicity (A = — 1) the replacement a O a’*' should 
be made in the Zeeman term, and a'^a — ^ Sp is con¬ 
served. 

Mean field theory. In the thermodynamic limit L —^ 00 
the bosonic state becomes macroscopically populated, 
and the mean field description (a) ~ \/nL exp(i(/)) be¬ 
comes exact [13) [l3- The 17(1) symmetry of the normal 
phase is spontaneously broken and the bosonic field ac¬ 
quires a nonzero, macroscopic mean value, implying su¬ 
perradiance. Note that integrating out the photon field in 
Eq. dU) leads to an effective long range electron-electron 
interaction of the form —g J^p k contrast to 

the short range interactions in low dimensional models. 


this long-ranged interaction facilitates the 17(1) symme¬ 
try breaking even at finite temperatures. 

The ground state is determined by either choosing the 
absolute minimum energy configuration!^ [l5| from all 
possible Q sectors, referred to as unconstrained solution, 
or by fixing the total number of excitations (25l - E3| in 
a constrained solution (i.e. grand canonical ensemble). 
The former situation is achievable in the presence of cou¬ 
pling to external noise or dissipation, inducing transition 
between various Q sectors, while the constrained solution 
corresponds to a situation in which polariton decay, aris¬ 
ing from coupling the cavity mode or the electronic exci¬ 
tations to other modes, is slow compared with the time 
required to reach thermal equilibrium at a fixed polari¬ 
ton number. The Q = 0 constraint is taken into account 
by adding — AQ to the Hamiltonian with A the Lagrange 
multiplier!^, . 

Within mean field theory, the fermionic excitation 
spectrum reads as 


Ea{p) 



( 2 ) 


where a = ± and, similarly to the case of classical light 0 
and to the one body picture, a gap A = opens in 
the spectrum due to the the quantized light - matter 
interaction. The mean field parameters, n, (j) and A are 
determined by minimizing the ground states energy 


L 


[uj — X)n 



( 3 ) 


where p = I/ttu is the ID density of states and IF is a 
high energy cutoff. The mean field equations, dE^ jdn = 
0 and OEq/ d\ = Q yield in the weak coupling limit 


\ P9 1 

oj — \ = -In 

2 




( 4 ) 


respectively. Both equations are to be solved in the con¬ 
strained case, while the unconstrained case requires only 
the first one in Eq. (|4]) with A = 0 and the second equa¬ 
tion should be ignored. 

The phase of the order parameter, (j) remains un¬ 
determined in both cases since the ground state en¬ 
ergy shows a Mexican hat in the space of Re (a) and 
Im(a). This degree of freedom corresponds to a gap¬ 
less Goldstone mode from the broken U{1) symmetry in 
the superradiant pha se of Eq. o, similar to the Tavis- 
Cummings model [ 23 . By tuning the phase (p adiabati- 
cally from 0 to 27r, one sweeps through the degenerate 
ground state manifold with no energy cost. The (p, (p) 
space can be mapped to the unit sphere with do, = 
aSf),{Sy)AS;)) = (Acos<).,Asin</.,up- A/2)/Eo(p). 
The winding number of the mapping defines the Chern 
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number, which reads for band a and given helicity of 
light A = ±1 as 

oo 27 r 

= J dp J • (5pda X = -a A. (5) 

— OO 0 


The Cn is analogous to the Chern number of the Floquet 


case 


ISj . The superradiant ground state manifold is thus 


topologically non-trivial, irrespective of the value of A. 

We emphasize again that the phase degree of free¬ 
dom from t/(l) symmetry breaking in a ID fermionic 
model arises due to the photon mediated infinitely long 
range interaction between electrons, while artificially cre¬ 
ated phases define similar Chern numbers in XY spin 
chain 2^ ■ A photocurrent can also generated along the 
edge as (j) = —eXl2Tr, depending on the value of A. 

Not only the bosonic field develops a macroscopic pop¬ 
ulation, but the time reversal symmetry is also broken, 
which is indicated by the finite electronic spin density 
perpendicular to the z direction as 


(S^) = -exp(=Fi'/>)- 


P_A^J2_W'^ 


A J 


( 6 ) 



FIG. 2. (Color online) The weak-coupling solution of the 
mean field equations, Eq. 0 is plotted for the constrained 
case, the low and high frequency asymptotes are indicates 
by dashed and dash-dotted lines, agreeing with the Floquet 
solution of the corresponding classical modelThe uncon¬ 
strained solution in Eq. © coincides with the constrained one 
in the high frequency limit. The inset shows the evolution of 
n for ujp = 1 as a function of the Zeeman coupling for the 
constrained and unconstrained cases (solid and dashed line, 
respectively). 


Unconstrained solution. The ground state is the ab¬ 
solute minimum e nergy c onfiguration over all possible Q 
sectors as in Refs. IJ,[ibj , since the total excitation num¬ 
ber is not fixed. The order parameter follows a typical 
weak-coupling, BCS-like form as 


n = 




to logarithmic accuracy. In the opposite, high frequency 
limit the solution coincides with Eq. 0. Due to the 
conserved number of total excitations, the superradiance 
of the cavity mode implies a finite photocurrent, which 
is also indicated by the finite value of A. The relevant 
energy scale of the problem separates the low and high 
frequency limits, which is signaled by a maximum in the 
photon occupation number, occurring at 


where the uj dependence is made explicit, and the gap, 
A = 2IFexp(—2 /p5 ) is independent of w. Since the 
electrons feel the photon mode through the opening of 
the gap, which is frequency independent, the physical 
response of the spin-Hall edge becomes also completely 
UJ independent. In contrast to the conventional Dicke 
model, a superradiant quantum phase transition occurs 
for arbitrarily small values of the coupling, g, while a 
critical coupling was required for the conventional, ho¬ 
mogeneous Dicke case. Eq. © indicates an infinite order 
quantum phase transition as a function of g as opposed 
to the second order transition of the conventional Dicke 
model[I^. Since A = 0 due to the lack of constraint, the 
photocurrent along the edge is zero. 

Constrained solution. Now we turn to the more rel¬ 
evant and interesting case with a fixed total number of 
excitations l25l427l | using both expression from Eq. (|4|). 


In the low frequency limit, the solution to Eq. o is 
obtained approximately as 






p'^guj 


\ P9^ J 


( 8 ) 


8 VF 

UJ = —^ exp 

pg 




(9) 


where the maximum of the photon number and current 
are Umax = P^gC/16 and jmax = -^eUmaxlT^P, as shown 
in Fig. [2l While the uj dependence features a maxi¬ 
mum for fixed g, the g dependence of these quantities is 
monotonic for fixed uj, shown in the inset of Fig. [2] The 
photocurrent behaves as 


(j) = J ^TT 2 y pguj J 
' 8eW^ / 4 A ^ 

- exp - for UJ ^ UJ. 

irpguj \ pgj 


( 10 ) 


The photocurrent arises from the constraint, which 
represents an effective magnetic field along the edge, 
inducing a photocurrent through the magnetoelectric 
effect [l|. In the low frequency limit, the current is 
pseudo-quantized as it grows proportional to uj (with 
log(I/w) corrections). This is similar to the Flo¬ 
quet case, where a perfectly quantized current with¬ 
out log-corrections follows from the adiabatic charge 
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pumping [^. The difference between the quantum 
and classical description arises because the condition 
of charge pumping that the Hamiltonian should change 
slowly and periodically in time is not satisfied here by 
the quantized photon field and quantum fluctuations do 
not allow for purely periodic time evolution and destroy 
quantization. This parallels to the breakdown of the 
classical equipartition theorem, which predicts vanish¬ 
ing energy per degree of freedom at absolute zero, which 
is violated by quantum fluctuations, which produce the 
zero-point energy of a quantum harmonic oscillator. 

In the high frequency limit, the current decays slowly 
as 1 /w with increasing frequency, which again agrees with 
the Floquet case. The ground state energy in this case is 
always above the absolute minimum energy state of the 
unconstrained case, although these approach each other 
fast in the lo ^ uj limit similarly to other quantities, 
because A vanishes as ~ — 1 /w with increasing frequency 
and A = 0 in the unconstrained case. By allowing for 
Q ^ 0 in the constrained case, the superradiant phase 
behaves similarly to Q = 0 and is accompanied by a 
finite photocurrent along the edge. For some uj and g 
pairs, the current may vanish, which indicates that the 
constrained solution in a given Q ^ 0 sector coincides 
with the absolute ground state. 

Due to the magnetoelectric phenomenon [l|, a finite 
magnetization develops in the z direction along the edge 
as (S^) = {j)lev, unlike in the unconstrained case. 

Photon spectral function. Having discussed the elec¬ 
tronic properties exhaustively, now we turn to the proper¬ 
ties of the photon mode in the superradiant phase, which 
contains the polariton spectrum and reveals the stability 
of the mean field solution. Its Green’s function is ob¬ 
tained by evaluating the Gaussian fluctuation correction 
of the order parameter around the mean-field value [i^. 
At T = 0, it reveals a zero energy Goldstone mode for any 
finite g due to the breaking of the U (1) svmmetrvj^. In 
the Nambu space, the Green’s function close to 

zero frequency becomes singular as 




n 


I a_p 
4 

gfp 

4 


a_p 

4 


-n 


a_p 

4 


( 11 ) 


and H is understood as H -|- iO'*'. The spectral function 

A(H) = Im Tr G{t^)/TT, measurable by the absorption 

coefficient of the cavity, features the Goldstone mode as 
2 

A(H —>• 0) = from the Gaussian fluctuations of 

the order parameter. The spectral weight of the Gold¬ 
stone mode vanishes at g = 0 due to the absence of phase 
transition there, and the spectral function assumes its 
non-interacting value as A(H > 0) = (5(0 — ui). The 
structure of the spectral function for vanishing O agrees 
with that of the finite-spin Tavis Cummings modelj^. 
In addition, a high energy mode is expected as a remnant 
of the bare cavity mode uj. This mode gets damped as 
^ TTg'^p for small A and is renormalized towards smaller 


frequencies before it hits the gap edge at O = 2A and 
merges with it. This optical gap is the amplitude mode 
(or Higgs like mode) where the spectral function displays 
a square-root singularity, as shown in Fig. [31 and this sets 
the threshold energy for continuum polariton excitations 
for fl > 2A. 



FIG. 3. (Color online) The contour plot of ln[a;A(fI)] reveals 
the polariton excitation energies and their spectral weights 
for IF/cJ = 100, Lop = 1. The white dashed line denotes the 
minimal excitation energy 2A, above which a continuum is 
formed. There is always a Goldstone mode at G = 0, omitted 
from these plots. 


Discussion. A single photon mode in Eq. ([T]) is real¬ 
ized in a quantum LG circuit 32| or is selected from a lad¬ 
der of cavity modes by placing a dispersive element into 
the cavity such as prism or nonlinear dielectric material 
(with wavevector dependent refractive index). The latter 
results in a non-equidistant ladder of excitations, which 
are eliminated by using dielectric mirrors (with high re¬ 
flectivity for the basic cavity mode, and high transmission 
for inharmonic frequencies). Tuning ui and g is possible 
by changing the dielectric constant or the position of mir¬ 
rors in the cavity, thus exploring various regions of the 
phase diagram. Multimode photon fields yield similar re¬ 
sults as was demonstrated in Refs. 14, 2^ This is 
most directly seen by integrating out the photon modes, 
which further enhances the effective electron-electron in¬ 
teraction. 

The SHI is realized using either condensed matter [i| 
or cold atomic [l^ [s^, [s^ settings. The electromagnetic 
field appears not only through the Zeeman term but via 
a vector potential. Neutral cold atoms in optical traps 
or chiral edge states (with spin parallel to momentum) 
[ssl do not couple to it, though. Otherwise, it can be 
incorporated into our calculations, and its main effect is 
the suppression of oj in Eq. 0, which is negligible for 
u! ^ vecy/pg/gespB-: where c is the speed of light, ^eff 
is the effective g-factor of the edge states and ps is the 
Bohr magneton. When a weak external field is coupled 
to the photon field, the phase of the order parameter can 
be tuned slowly and the quantized Ghern number can 
be probed through a Thouless typejs^ adiabatic charge 
pumping. 

The superradiant phase transition occurs for arbitrary 
small coupling between the photon field and the spin- 
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Hall edge, and is characterized by a Chern number ±1 
in the space of the ID Brillouin zone and the phase of 
the order parameter. The unconstrained solution, cor¬ 
responding to the absolute ground state of the system, 
does not carry a net electric current. By keeping the to¬ 
tal number of excitations fixed, a finite electric current 
carrying state develops, which behaves analogously to its 
Floquet counterpart. This remotely parallels to the idea 
of quantum time crystals [s^, which do not occupy the 
absolute ground state of a system but appear in an ex¬ 
cited sector 37| . The conventional Dicke model exhibits 
quantum chaos 1^ away from the thermodynamic limit 
as a precursor to quantum phase transition at L —>■ oo. 
It would be interesting to investigate whether and how 
quantum chaotic behaviour sets in in the present model 
for finite L. Extending our work to higher dimensional 
systems, which are suspected to become Floquet topolog¬ 
ical insulators such as graphene[38| or HgTe/CdTe quan¬ 
tum wellsQ or their cold atomic realizations, promises 
to be a fruitful enterprise. 
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